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Abstract 

In the history of mechanics, there have been two points of view for 
studying mechanical systems: The Newtonian and the Cartesian. 

According the Descartes point of view, the motion of mechanical 
systems is described by the first-order differential equations in the N 
dimensional configuration space Q. 

In this paper we develop the Cartesian approach for mechanical 
systems with three degrees of freedom and with constraint which are 
linear with respect to velocity. The obtained results we apply to dis- 
cuss the integrability of the geodesic flows on the surface in the three 
dimensional Euclidian space and to analyze the integrability of a heavy 
rigid body in the Suslov and the Veselov cases . 



1 Introduction. 



In " Philosophiae Naturalis Principia Mathematica" (1687), Newton considers 
that movements of celestial bodies can be described by differential equations 
of the second order. To determine their trajectory, it is necessary to give 
the initial position and velocity. To reduce the equations of motion to the 
investigation of a dynamics system it is necessary to double the dimension 
of the position space and to introduce the auxiliary phase space. 

Descartes in 1644 proposed that the behavior of the celestial bodies be 
studied from another point of view. These ideas were stated in "Principia 
Philosophiae" (1644) and in "Discours de la metode" (1637). According to 
Descarte the understanding of cosmology starts from acceptance of the initial 
chaos, whose moving elements are ordered according to certain fixed laws and 
form the Cosmo. He consider that the Universe is filled with a tenuous fluid 
matter (ether), which is constantly in a vortex motion. This motion moves 
the largest particle of matter of the vortex axis, and they subsequently form 
planets. Then, according to what Descartes wrote in his "Treatise on Light", 
"the material of the Heaven must be rotate the planets not only about the 
Sun but also about their own centers... and this will hence form several small 
Heavens rotating in the same direction as the great Heaven." [2]. Thus the 
equation of motion in the Descartes theory must be of the first order 



Hence, to determine the trajectory from Descartes's point of view it is nec- 
essary to give only the initial position. 

In the modern scientific literature the study of the Descarte ideas we can 
find in the monographic of V.V. Kozlov [2] in which the author give the 
following result. 

Theorem 1.1 The manifold y = u{x,t), where u is a covector on Q is an 
invariant manifold for the canonical Hamiltonian equations with the Hamil- 
tonian H{x,y,t) if and only if field u satisfies the Lamb equation 





dtu{x,t) + {rotu{x,t))v{x,t) 





where (rotu) = d^u — {d^u)^ ia skew- symmetric n x n matrix, 



v{x,t) = dyH{x,y)\y=uix,t), h{x,t) = H{x,y,t)\y=u{x,t) 



(1.3) 
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From the physical standpoint, equations (II .ip . fll.2p and fll.3l) describe the 
motion of the colhsionless medium: particles moving along different trajec- 
tories do not interact. 

In [2] affirm that "solving dynamics problem is possible inside the con- 
figuration space". For this it is necessary to solve Lamb equations which 
is a system of partial differential equations on Q, and then, using (11.31) to 
calculate the vector field v from the solution of the Lamb equation to solve 

(HH). 

In [3j we developed the Cartesian approach for mechanical system with 
configuration space Q and with constraints hnear with respects to velocity. 
The aim of the present paper is to develop the results obtained for mechanical 
system with three degrees of freedom in the particular case in which Q = 

is the three dimensional Euclidean space and Q = .SO (3) is the special 
orthogonal group of rotations of 



2 Cartesian vector field on three dimensional 
Euclidean space 

Let ¥? be the three dimensional Euclidian space with cartesian coordinates 

X = (Xi,X2,X3). 

We consider a particle with Lagrangian function 

l = Im'-u{x) 

and constraints 

(x,a) = (2.1) 

where (, ) denotes the scalar product in K^, x = (xi, ±2, X3) and a(x) = 
(ai(x) , a2{x), 03 (x)) is a smooth vector field in 

It is well known that the equations of motions can be deduced from the 
d'Alembert-Lagrange principle [9] 

x= [/'x + /ia(x), . 
(a,i) = 0,, 

where /i is the Lagrangian multiplier, f/x = {Ux^, Ux2, Ux-j), Ux^ = dx^U. 
In [3] we introduce the following definition 
Definition 1 
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We say the smooth vector field v{x) = {vi{x), V2{x), fa (a;)) is the Carte- 
sian vector field for a constrained particle in with the constraints \2.1\] 

[w{x) X rotv(x)] = A(x)a(a:). (2.3) 
where [ x ] denotes the vector product in E"^, 

and A is a function: 

A = ([v(x) X rotv(x)], a(x)). 
a r 



By a simple computation from (12.31) we can see that 

(a(a;), v(x)) = 0, 
(a(a;), rotv(x)) = 0. 



(2.4) 



Corollary 2.1 Let v the Cartesian vector field. Then the following relations 
hold 



\\Hx)\\'^)^ + K{x)a{x) ^2.5) 



^2 



[x, a) = 0, 

The proof it is easy to obtain in view of the equality 

X = (^||v(x)|n^ + [v(x) X rotY{x)] 

which is deduced after derivation the differential equations generated by the 
vector field v 



i = v(x) (2.6) 

The system (12. 5 p can be obtained from the Lagrangian equations with La- 
grangian function 




where v is a Cartesian vector field. 
Definition 2 

The study of the behavior of the constrained particle in E^ by using the 
equations h2.2^) or h2. 6\) . I[27^) or l{2.5\) say the Classical, Cartesian and La- 
grangian approach respectively. 
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We illustrate the above ideas in the following example 
A non-holonomically constrained particle in R'^. 



Consider a particle with the kinetic energy T = || |x| p and non-holonomic 
constraints 

xi + a{x^)x2 = 

This instructive academic example, in the particular case when ^(xs) = x^ 
due to Rosenberg [S]. This example was also used to illustrate the theory in 
Bates and Sniatycki [1]. 

The Descartes approach in this case produces the vector field v : 

V = X2{a{x3)dxi - dx^) - \3dcc3 
and condition fl2.4l) for this case takes the form 

(rotv, a(x)) = ^d^aii^ + 0-^)^2) + {a{x3)d^^X3 - d^^X3)\2 = 0. 

We shall study the case when this relation holds in view of the equalities 

A 



A, 



A3 = b2{x3), 



for A an arbitrary constant and 62 an arbitrary function on X3. 

The equations generated by the vector field v in this case can be written 

as 

a{x3)A 



Xi 



X2 



A 



(2.7) 



^l + a^ixs) 

\^ X3 = -b2{x3) 

The all trajectories of these equations are easy to obtain. 

The Lagrangian approach produces the following differential equations 



Xi 



-b{x3)d^. 



Aa{x3) 



X2 = 6(2:3)5^.3 ( 



a/1 + a^ix-s) 
A 



a/1 + a^ixs) 



^3 = d^3^b^{x3) 
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Corollary 2.2 All the trajectories of the equation of motion of the con- 
strained Lagrangian system 

< L = - U{x3), {xi + d{x3)x2 = 0} > 

can be obtained from ^2. 7| j with b{x3) = ±a//i + U{x3). 

In fact, the equations of motion obtained from the D'Alembert-Lagrange 
Principle are 

' xi = 

X2 = a(x3)/i 
^3 = dx3U{x3) 
Xi + d{x3)x2 = 



Therefore, 



hence. 



d ^, ... dd[x3) . . 
— (xa - a(x3)a;i) = — XgXi 



A 

X2 



where A is an arbitrary constant. 

On the other hand from the equation 

X3 = dx.^U{x3) 

we easily obtain X3 = ^^/2{U{x3) + h), where h is an arbitrary constant. 

Finally by considering the constraints we deduce the system of the first 
order ordinary differential equations (12. 7p . In this example the Descartes, 
the lagrangian and Classical approach coincide . 

Below we determine the Cartesian vector field for a particle on the surface 
in E^. 

First we introduce the ve ctor fields X, ,Y, Z which are characteristic 
elements of the 1-form 

fl = ai{x)dxi + a2{x)dx2 + a3{x)dx3 

X = a3dy - a2(9^ 

Y = aid, - a3dx (2.8) 
Z = a2dx — CLidy 
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Clearly, the most general element of the given 1-form Q is 

V = wiX + W2Y + W3Z 

hence 

v(x) = [a(x) X w(x)] (2.9) 

where w(x) = {wi{x), W2{x), 103 (x)) is an arbitrary smooth vector field 
which we shall determine in such a way that (12. 4p takes place. 
By using the identity 

rot[a(x) X h{x)] = [a, b] + divh{x)a{x) — divSL{x)h{x) 

where [a, b] is the Lie bracket of the smooth vector field a and b, one can 
prove the following assertion 

Corollary 2.3 The condition ^K^,^MIM) ca?^ be written as follows 

div{[a{x) X [a{x) x t«(x)]]) = {[a{x) x rota{x)], w{x)) (2.10) 
Proposition 2.1 Let us suppose that the vector field a is such that 

[a{x) X rota{x)] = (2.11) 
then the Cartesian vector field exist i and only if 

a{x) = fx{x) (2.12) 
for a certain smooth function f. 
Proof 

From (I2A0D . fl2TTD follows that 

[a(z) X [a(x) X w(a:)]] 
is a solenoidal vector field, hence 

[a(x) X [a(a:;) x w(a;)]] = rot'W{x) 
for arbitrary vector field W, thus the following representation holds 

a(x),w(x)) rotW(x) 



w = — ,, , ,,,0 — a x + 



lalxllp lla(x)|p 



as a consequence 

(a(x),rotW(x)) = 0. (2.13) 

Clearly if (a(x), rotW(x)) 7^ then the Cartesian vector field does not 
exist if (12. lip holds. If we choose 

a(x) = /x(a;), W(x) = ^Gxl/, $) 

when $, G are an arbitrary smooth functions, then we obtain that (12.131) 
holds identically and a consequence the vector w takes the form 

Corollary 2.4 The Cartesian vector field for a particle in ¥? which is con- 
strained to move on the surface 

f{x) = c,c^O (2.15) 

generated the following differential system 

x=uix)i\\fxix)\\^^xix) - ifxix), ^xix))fxix)) (2.16) 

Corollary 2.5 The Lagrangian approach for a particle in which is con- 
strained to move on the surface 1^2. 15\} produces the following differential 
equations 



x= ^{li^\x)\\fx{x)\mx{x) X <!>x{x)]\\') +A{x)fx{x) (2.17) 



Corollary 2.6 // there exist a function G and $ such that 



1 1 [fx{x) X $^(x)] IP = ^ ^{f, $) g (2.18) 



Then the equations \2.11\) take the form 



i = Ao(x)/x(x), \{x) = hf{f) + K{x). (2.19) 
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If one introduce the matrix A{x) : 

(fxixi fxiX2 fxiX3 
fx\X2 fx2X2 fx2Xz I (2.20) 
fxiX3 fx2X3 fxgX^-, 

then one checks, that the equations (I2.17P may be written as 

iA{x)v{x),v{x)) ^ 

where v is the Cartesian vector field generated the differential equations 
( I2.16p . The differential equations (12.191) determined the geodesic flows on 
the surface (12.151) and admits the energy integral 

M? = 2h{f). 

If there is an additional flrst integral, functionally independent with the 
energy integral , then the geodesic flow is integrable. 

In order to study the integrability of the geodesic flow on the given surface 
we introduce the following functions which we determine from (I2.16P 



ll/x||||[$xXx ]||V 
(/x, $x) 

||$X||||[XX X]||\2 



ll[xx<I.x]|| 
In view of (12.18^ . it is easy to show that 

F, = 2h{f), J = 1,2,3. 



if (/x, $x) 7^ 0; 



if (/x, $x) =0, $x 7^ «:(x)x. 



3 Integrability of the geodesic flow on the ho- 
mogeneous surface. 

We now consider the surface 

/ fix) = c, c ^ 0, , . 

\ (x,/x(x))=m/(x). ^^-'^ 
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which wc will call the homogeneous surface of degree m. 

From the Euler formula follow that c = is the unique critical value of 
/. hence for c 7^ the function 

^ = ||/x(x)|r>0 

on the given homogeneous surface. 

Taking (3.1) into account we deduce the relations 

A(x)x^ = (m-l)/^(x), (3.2) 
(x,5x(x)) = 2{m-l)g{x) (3.3) 

Below we use the following notation 









F 


{F, G, H) = 






Gxz 






Hx2 





Clearly, if F, G, H are independent functions then {F. G, H} ^ 0. 
The integrability of the geodesic flow on the homogeneous surface we shall 
study in the following two cases 

{f,9,r'} = 0, (3.4) 
{f,9,r'} + (3.5) 

where r"^ — x\ ^ x\ ^ x\. 

We analyze the first case. We study only the particular subcase when the 
homogeneous surface is such that 

ll/x(a;)|P=5(/,r). (3.6) 

Hence, in view of (3.3) we give 

mfdfg + rdrg = 2{m-l)g{f,r). (3.7) 

We assume that the arbitrary function $ is such that 
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thus the differential equation generated by the Cartesian vector field and 
second order differential equations of the geodesic flows under the indicated 
condition take the form respectively 

X = z/(x)(5(x- m//x), (3.8) 

^^rna^ (3.9) 

where 

iy'g{f,r){g{f,r)r'-m^f) (3.10) 

Proposition 3.1 The geodesic flow on the homogeneous surface under the 
assumption ^3.6\) is integrable 

Proof 

First we observe that there is the function ly such that fl3.10p holds , i.e., 



2hif)g{f,r) 
g{f, r)r'^ — rn^p ' 



hence exist the additional first integral F2 which in this case takes the 
form 

F2 = aU.mUxUm' ^ g(^f^r)\\[Ux) X x]|p = 2m'fh{f). 

The particular class of the study homogeneous surface are the following. 
If m = 1 then (x, ^fx) = 0, in particular this relation holds if 

g = '^{-) 

r 

A concrete example we obtain from the celestial mechanics |5]: 

/(x) = r+(b,x) = c, c^O, (3.11) 
where b = (fei, 62, 63) is a constant vector field. In this case we have 

r 
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The first integral for tliis particular case are 



'2^ + ||b|r-l)||[xxi]||^ = 2A(/) ^^-^^^ 



r 

It is interesting to deduce the equations of motion of a particle constrained 
to move in the surface (13. lip with the subsidiary condition that there is a 
nonzero constant vector field c = (ci, C2, C3) : 

(x,c) = 0, (b,c) = ^(/x(x),c) = 0, (3.13) 

from the Lagrangian approach. 

By choosing the function $ as follows 

$ = (C3 - C2)xi + (Ci - C3)X2 + (C2 - Ci)x3 

and introducing the new time a as 

da = ( h Oi + 02 + 03)"^ 

r 

and letting the prime denote differentiation with respect to a, we have that 
the equation generated by the Cartesian vector field can be written as 

= n^lM^) X (3.14) 



By considering that 

rot[fx X c] = - 
r 

we obtain that the Lagrangian approach generated the second order differ- 
ential equations 

/(x)x 



igi 12^3 



Thus , if 

then we obtain the well known equations 



//X , . 

X =-— (3.15) 



12 



These equations admit the following first integrals 



,2 2 
= - + 
r 



llb 


2 


- 1 




C 


2 



[x' X c] 



;-+b) 

r 



(3.16) 



[x X x'] = c 

These relations are easy to obtain from (13.181) . 

The equations (I3.18p . after the orthogonal transformation 



(b,x) 



take the form 



V 



([c X b],x) 
||c|| ||b|| 

V 



c 



(c,x) 



(3.17) 



I C = o. 

These equations describe the behavior of the particle with Lagrangian func- 
tion 

constrained to move on the one curve of the family of conies 

f{x) = ^/e + v' + \\b\\^=\\c\\' 

The differential equations generated by the Cartesian vector field in this co- 
ordinates can be represented in Hamiltonian form with Hamiltonian function 

fm 

We now turn to the study the particular case of the homogeneous surface 
with the condition g = g{f,r). 
If 
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then after the change 



we obtain 



Finally, if 



2(m.-l) , f , 



then after the change 
we deduce the equation 



which show that this case is equivalent to the first case study above. 

We have aheady studied the case in which {/, f?,?"^} = 0. Now we begin 
to study the case in which the functions /, are independent. Hence 

{Lg,r'}^0. (3.18) 
Under this assumption we obtain that 

Xj ^ Xj{f,g,r^), j = 1,2,3 

thus we deduced that 

Proposition 3.2 // there exists the functions $ and G such that 

<^ = g,r'), G = G(/,$) 

gu = Gf{f,^) (3.20) 
Wfximxx fx\\'i^'^2h{f) 

then the geodesic flow on the homogeneous surface of degree m > 1 is inte- 
grable. 
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Proof 

We prove this assertion only for the case when 

($x, /x) = 0, (3.21) 

thus the surface $ = ci is orthogonal to the given homogeneous surface. 
Under this assumption we obtain that the differential equations deduced 
from the Cartesian and Lagrangian approach are respectively 



X = z/(x)$x(x) 
X = Ao(x)/x(a:) 



(3.22) 



where An can be determined as follows 



|2 



Ao(x) = /^^(/x,5x(^||$.|P)) = hfif) - ^i^(/x,j.x) 
After derivation the function 

|<l>xllllix X xl 



||[xx$x] 

along the solutions of (13.221) we obtain 



dt 



($x,5x(/^1l$x|r)). 



which is equal to zero in view of (13.201) . (13.211) thus the function F3 is the first 
integral of the geodesic flow. 

Clearly, in order to assess the integrability of the geodesic flow in this 
case we need first to check whether function u exists such that (13.201) holds. 



4 The geodesic flow on the quadrics and the 
third-order surface 

In order to illustrate the above ideas we consider the algebraic surface of 
degree three: 

f{x)=XiX2X3 = C, c 7^ 0. (4.1) 
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This case was examined already by Riemann in his study of motion of a 
homogeneous hquid elhpsoid. More exactly, Riemann examined the integra- 
bility of the geodesic flow on (14.11) . 

In 1^ the author state the following problem. 

" Is it true that the geodesic flow on a generic third-order algebraic surface 
is not integrable?. In particular I do not know a rigorous proof of non- 
integrability for the surface (14.11) " 

By considering that in this case 

g = {xiX2f + {xixsf + (^3X2)^ 

thus the functions f,g and are independent. The dependence Xj = 
Xj{f, g, r^), j = 1, 2, 3 we obtain as follows. 
We introduce the cubic polynomial in z : 

P{z) =z'- r'z' + gz~f = {z- xl){z - xl){z - xl), 

and by using Cardano's formula we obtain the require dependence. 

In order to construct the Cartesian approach in this case first we observe 
that the surface 

^{C,v,0 = ci 

where 



1 1 1 

? ~ 2^^1 ^2/' V — 2^^3 ^iJ' s — 2^'^2 ^3) 

is orthogonal to surface (14. ip . Thus the differential equations generated by 
the cartesian vector field are 



Xl = i^($5 — ^n) Xl 

i;2 = zy(<l>^ - $5) X2 (4.2) 
X3 = H% - ^() ^3 

To determine the existence the solution of (I3.20p or, what is the same, 

2 _ 2/i(/) ^ 
11^x11 G'^(/,$) 

is for us an open problem. 

Now we study the subcase when the given surface is such that 

fix) = ^{bixl + b2xl + hxl) (4.3) 
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First we state and solve the following problem. 
Problem 1 

Let X, y, Z are the vector fields ([2l]), flTOD . 

We require to determine the function / in such a way that these vector 
field formed a three dimensional Lie algebra. 

The solution of this problem it is easy to obtain in view of the equality 

Ti = A(x)T2, (4.4) 
where A is the matrix given by the formula ([2]) and 

Ti = col{[Y,Zl [Z,Xl T2 = col{X, Y, Z) 

and by using the Bianchi representation 

T3 = B{x) T4, (4.5) 

where T3 = col(\U, V], [V, W], [W, U]), T4 = col{U, V, W), where U, V, W 
are the vector fields, B is the matrix: 

/ a h \ 
B{x) = 61 

V 62 -a / 

and a, bi, 62, 63 are certain constants 

From fl4.4p and 04.51) after integration we obtain the class of functions 
which generated the three dimensional Lie algebra. 

1. f = hx' + hy' + hz' 

2. f = b,x^ + a{y^-z^)+2byz 

3. f = 2byx + bsz^ 

4. f = by'^ + 2bizx 

We construct the Cartesian vector field for the first case. 
In view of the relation 

g = bjxf + bjxl + 63X3 

we observe that /, g, are independent functions. The following equalities 
it is easy to obtain: 
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-2{b2 + b3)f + g 








-2{b^ + bs)f + g 




-b,){b2-bs) 


b2bir^ 


-2(62 + 61)/ + ^ 



Notice that 
where 



{b3-b2){bs-b,) 



™>^3 

3 2 /■ 1 



(4.7) 



^2 -^1 -^3 

is an orthogonal surface to the given surface we obtain that the differential 
equations generated by the Cartesian vector field in this case can be written 
as 



Xi 



Xl 



X2 = —{%vbi - ^cC^s) 

X2 



Xl = — ($^e62-$cC^i) 
X3 



(4.8) 



Now we introduce the elliptic coordinates in : 

(Ai + 6r')(A2 + 6r')(A3 + 6r') 



X 



X 



{bT\- b^')ib^] - b^') 

, {X, + b^'){X2 + b^'){Xs + b^') 

(62-^-6r^)(62-^-63-i) 
^_(X, + b^'){X2 + b^'){Xs + b^') 



(4.9) 



ib3'-b2'm'-bi') 
where Ai, A2, A3 are the roots of the cubic polynomial in w : 

-w' + {xl + xl + xl - bT' - 62-' - b^')w' + ixjib^' + 63 ') + + b^') 

+:r2(6-i + 6-1) - b^^b^^ - b:[%^ - b^%^)w 
+b^%%\bixl + b2xl + b^xl - 1) = 
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In this coordinates we obtain 



\x\ 



where 



9u 



922 



933 



(Ai-A2)(Ai-A3) 

4(Ai-6r')(Ai-62')Ai-63'^ 
(A2-Ai)(A2-A3) 



4(A2-6r')(A2-62^)A2-63 . 

(A3 — A2)(A3 — A3) 



(4.10) 



4(A3-6r^)(A3-62"')A3-fc3 
The differential equations ( 14. Sp in elhptic coordinates can be transformed to 
the form 

r w = o 



2u 



u 



{{h2 - ?'3)$cC + {bi - h2)^r,V + % - = ^i(Ai, A2) 



616263 

T^ir- + T- + - h)hi^cC + (61 - h2)h^^r] + (63 - &i)&2$^0 = ^2(Ai, A2) 

61 62 63 Oi 62 63 
^ = Ai + A2 + A3, V = A1A2 + A1A3 + A2A3, w = A1A2A3, 

and by putting 

A3 = 0, 

after some calculations we deduce the planar system 

^i(Ai, A2)Ai - ^^2(Ai,A2)) 



Ai 



A, 



^i(Ai,A2)^A2 -^^2(Ai, A2)) 
Ai — A2 



In order to deduce the differential equations for the geodesic flow by using 
the Lagrangian approach first, it is necessary in the first place obtain the 
solution of the equations (13.201) . 

The integrability of the geodesic flow on the quadric (m=2) by using the 
classical approach, was proved by Jacobi and Chasles. 
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5 The geometrical and physical meaning of 
the Cartesian vector field 



The purpose of this section is to determine the geometrical and physical 
meaning of the Cartesian vector field constructed above. 

Hertz's principle of least curvature is a special case of Gauss' principle, 
restricted by the two conditions that there be no applied forces and that all 
masses are identical. (Without loss of generality, the masses may be set equal 
to one.) Under these conditions. Gauss' minimized quantity can be written 



N 



cPx^ 



2 



The kinetic energy 



2 



is also conserved under these conditions 

Since the line element ds^ in the 3A/"-dimensional space of the coordinates 
is defined 

ds^ = 2Tdf ^^ = 2T 
dt^ 

by considering the conservation of energy we obtain 

ds^ 

Dividing Z by 2T yields another minimal quantity 



N 



I ^ '^'^ 1 2 



Since v K is the local curvature of the trajectory in the 3 A^- dimensional space 
of the coordinates, minimization of K is equivalent to finding the trajectory of 
least curvature (a geodesic) that is consistent with the constraints. Hertz's 
principle is also a special case of Jacobi's formulation of the least-action 
principle. Curvature refers to a number of loosely related concepts in different 
areas of geometry. In mathematics, a geodesic is a generalization of the 
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notion of a straight line to curved spaces. Definition of geodesic depends on 

the type of curved space. If the space carries a natural metric then geodesies 
are defined to be (locally) the shortest path between points on the space. 

Below we restricted to the case when the configuration space is the three 
dimensional Euclidean space with Cartesian coordinates x = {xi, X2, X3). 
The geodesic flow on the surface f{x) — c is determined by the second- order 
differential equations 



cPx^ _ _ (A(a;)x, x) . _ . ^ „ 

^^2 -^Z-.' ||j^||2 J--L>A>i 



which, by considering the energy integral, can be written as follows 



d^x^ ~ r ~ {A{x)t, t) . 



where 



Clearly that 



as 
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ll/xl 

The Hertz's Principle of Least Curvature and problem on the determi- 
nation the principal directions on the surface lead us to state the following 
problem. 

Problem 2 

Determine the 

under the conditions 



extre'mum{A{x)T, r) 



||r||2-l =0 
, (/x,t) = 

Solution 

Note that in this case the Lagrangian function is 

L^{A{x)T,r)+a{f^,T) + zi\\T\\'-l) 
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where a and z are the Lagrangian multipher and computer 



r dh 

dL 



:0, j = l,2,3«(A(x) + z/)T^ + a/^ = 
0^(/x,r) = 



dL 

. dz 



^ = ^ llrl 







where = coI{ti, T2, ts) and / is the diagonal matrix: I = diag{l, 1, 1), 
from the first group of equations, we deduced the following equahties 



and 
if 



T^{x) = -aiA{x)+zir'f^{x)) 



(5.1) 
(5.2) 



det {A{x) + zl) ^ 0, 
where x = col[ri, r2, r^, a) and is the following family of matrixes 



In view of that the vector x is non-zero vector then from 05.21) one can 
deduce that 



^ fxixi ~\~ Z 


fxiX2 


fxiX3 


fxi 


\ 


fxiX2 


fx2X2 "1" ^ 


fx2X3 


fx2 




fxixs 


fx2X3 


fx3X3 + Z 


fx3 




\ fxi 


fx2 


fx3 





/ 



det R^ = 



(5.3) 



In order to establish the relation between the vector field with components 
given by (15.11) and Cartesian vector field (I2.16P we introduce the family of 
vector fields v, : 



= {{A{x) + zl) Vx,4), = {d^,, 9^2, a 



X3J 



(5.4) 



where ^ is a complex parameter. After some calculations one can prove 
that Vz admits the representations 
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fxixi ~l" Z 


fxiX2 


fx 1X3 


fxi 


fxiX2 


fx2X2 ~l~ ^ 


fx2X3 


fx2 


fxix:i 


fx2X'i 


fx^xs ~l" ^ 


fx'i 




9x2 


9x3 






det {A + zl) 



(5.5) 



which are equivalents equivalent to (15 .4^ . where A/ = 9x^xif +9x2X2 f +9x^x-j,f ■ 
and Xi, X2, X3 denote the vector fields: 



^1 — {fx, 9x) 

X2 = {gx,9x) - Af{f^,dx), g=\\fx\\'^ 
X3 = dety4(2;) v^\^=o 

We now introduce the function 

F(.)^./(..)^^™ + « + ^^™ 



(5.6) 



(5.7) 



det {A + zl) 
which in view of (15. 5p may be written as 

F{z) = det (A + zl) det 
Corollary 5.1 If the function f is a homogeneous function of degree m then 



X3 = g\x)K{x){x,dx), g=\\fx 



(5.8) 



where K is the Gaussian curvature of the homogeneous surface which one 
can calculate as follows 



K{x) 



fifxi, fx21 fxs} \{ -L \- 

(m — l)g'^ ' 

f{fl fx21 fx^} 



(5.9) 



if m = 1; 



xig^ 



Below we shall study only the case when the function / 
idfiX2)y-AdfiX,)dfiX,)>0 



(5.10) 



Clearly, under this assumption the function F has two different real roots 
which we shall denote by zi, Z2 
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Corollary 5.2 Let Vi, V2 are the vector fields such that 
then the solution of the problem 2 are the vector fields 





J = 1,2 



(rW, r(2)) =0 



The proof it is easy to obtain. 

Proposition 5.1 Let f{x) = c, c ^ be the homogeneous surface which 
satisfies l\5.10\) . 

Then the most general vector field tangent to the given surface admits the 
development 



One can check direct from the above that the most general vector field 
tangent to the given homogeneous surface can be written as 



(5.11) 



where fii, /i2 are arbitrary smooth functions. 



Proof 



v(x) = ai{x)v^^\x) + a2{x)v^'^\x) 



where ai, 02 are arbitrary smooth functions. 
A brief calculation show that 



v{x) 



/x(a::) + \2{x) gy^{x) + \3{x) x 



(5.12) 



where 



Ai(a;) 
A2(a;) 
A3 (a;) 

Uj (x) 



ai{x){zl - Af zi) + a2{x){z1 - Af Z2) 

ai{x) zi + 62 (x) Z2 

g'^{x)K{ai{x) + a2{x)) 

hj det {A{x) + ZjL), j = 1,2 



One can see that the equation 
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^7Ai + (/x, <?x)A2 + m/A3 = (5.13) 

holds identically. 

To complete the proof, we show the equivalence of (15. lip and (I5.12p . 
(I5.13p . Indeed, using (15.131) we obtain that 

A,--( - ), 

inserting into (I5.12p and introducing the notations 

/ N A3 . X A2 
^i[x) = — , jj,2{x) = — 

9 9 

we get flSTTT]) . 

Proposition 5.2 The vector field Ii5.11\) is Cartesian vector field if the fol- 
lowing relation holds 

{f,g,r}{drfii-dgfi2-j) = (5.14) 

From the definition we obtain that the given vector field is Cartesian (see 
definition 1) if the following equality takes place 



(/x,(rot([/x(x) X [/x(x) X (^i(x)^x(2:)+/i2(x)x))]]) =0 (5.15) 
which is equivalent to ( I5.14p . 

From (I5.14P after straightforward calculations we can prove the following 
assertion 



Corollary 5.3 Let us suppose that 1^5. 14\ ) holds, then the vector field \5.12^) 
admits the representation 



v(x) = [/x X [/x X <l>x]] (5.16) 

where k and $ functions such that 

r «:(/,$) if{/,^,r}^0; 

K=\ (5.17) 
[ /^i 7" + if {f,9^'r} = 0; 
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and 

r <l>x(/,^,r) if {/,^7,r}^0; 
= <^ (5.18) 

[x if {/,^7,r} = 0; 

By comparing fl5.16l) . fl5.17p with (12.161) we obtain that the solution of 
the problem 2 under the condition (I5.17p . (15.181) coincide with the Cartesian 
vector field constructed in the section 3. In such a way we obtain the physical 
and geometrical meaning of the constructed above Cartesian vector field. 

Remark. 

The physical and geometrical meaning of the Cartesian vector field for 
the case when the given vector field a : 



(a,rota)^0 (5.19) 

can be obtained analogously to the case study above by considering that 
under condition (15.191) the equation deduced from the Lagrangian approach, 
can be written as follows 

(yl(x)x, x) 



-a 



where 



ditti ^{dia2 + d2ai) ^{dias + d3ai) 

A{x) = ( ^{dia2 + d2ai) ^202 |(c^2a3 + <93«2) 

lidia^ + dsui) 1(^203 + 5302) 9303 

The problem 2 in this case we can state as follows 
Problem 3 
Determine the 

extremum{A{x)T. r) 

under the conditions 

' Ikll'-i = 

(a,r) = 
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6 Descartes approach for non-holonomic sys- 
tem with three degree of freedom and one 
constraints . 

Our goal in this section is to extend the Cartesian approach developed above 
for natural mechanical system with configuration space 

Q, dim Q = 3 

in this space the metric (kinetic energy) 



1 ^ 

2^ = o E GAx)x^i' 



_ ikiXjX'X = - \\x\\ 

j,k=i 

allows calculating the rot of the vector field on Q. The invariant definition 
of rotv we can find in [2]. If we assign a covector field p = {pi, P2, Ps) with 
components 

3 

Pj = Y.G,,{x)v\x) (6.1) 

k=l 

to the vector field v = {v^{x), f^(x), v^{x)), then the components of rotv we 
can write explicitly 

111 

rotv{x) = {—={d2Pz - d3P2), -j={d3Pi - diPs), —={diP2 - d2Pi)) 
vG V Ct V Ct 



where G = det {Gkj{x)). 

The vector field fl2.9p for this mechanical system we shall represented as 
follows 

X = [a(x) X (Ai h{x) + A2c(x))] (6.2) 
where a, b, c are the independent smooth vector in Q, i.e., 

T = (a, [b X c]) ^ (6.3) 

and Ai, A2 are smooth function which we determine as a solution of the 
equation 

(a,rot[ax (Aib + Aac)]) =0 (6.4) 
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The Lagrangian approach produces the following second-order differential 
equations 

We shall illustrate this case for the Chaphguin-Caratheodory sleigh and 
for the heavy rigid body in the Suslov case. 
The Chapliguin-Caratheodory sleigh 

We shall now analyze one of the most classical nonholonomic systems : 
Chapliguin-Carathodory's sleigh [15J. The idealized sleigh is a body that has 
three points of contact with the plane. Two of them slide freely but the third, 
A, behaves like a knife edge subjected to a constraining force M which does 
not allow transversal velocity. More precisely, let yoz be an inertial frame 
and ^ Ar] a frame moving with the sleigh. Take as generalized coordinates 
the Descartes coordinates of the center of mass C of the sleigh and the angle 
X between the y and the ^ axis. The reaction force M against the runners 
is exerted laterally at the point of application A in such a way that the 77 
component of the velocity is zero. Hence, one has the constrained system M 
with the configuration space X = x M^, with the kinetic energy 

and with the constraint 

ex + sin xy — cos xz = 0, 

where m is the mass of the system and Jc is the moment of inertia about a 
vertical axis through C and e = \AC\. Observe that the "javelin" (or arrow 
or Chapliguin's skate) is a particular case of this mechanical system and can 
be obtained when e = 

To apply the Descartes approach for this system, first we determine the 
vector b and a in such a way that the determinant T 7^ 0. In this subcase, 
we achieve this condition if 

a = (e, sinx, — cosx) h = {0, cos x, sinx), c = (1,0,0). 

Under these restrictions we obtain that T = 1 and it is easy to show that 
the vector field v takes the form: 

V = X-ii^Bx + esinxBy + ecosxBz) — X2{cosxBy — sinxB^). 
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The Descartes approach produce the differential equations 

X = X3{x,y,z,e) 

y = X2{x,y, z, e) cosx — eX^sinx (6.6) 
z = A2(x, y, z, e) sin x + eX^ cos x 

where A2, A3 are solutions of the partial differential equations 

sin x^JdzXs + emdyX2) + cos x^JdyX-s — tmdzX2) — m{dxX2 — eAs) = (6.7) 

where J = Jc + e^m. 
Clearly, 

I |v| 1^ = {Jc + me'^)Xl{x, y,z,e) + mXl{x, y,z,e) 

Hence, for the arrow (e = 0) we have 

X = X3{x,y,z,0) 
y = X2{x, y, z, 0) cosx 
i = X2{x, y, z, 0) sin a; 

Jc(sinx9^A3 + cosxdyXs) — mdxX2 = 
Clearly, the equation (16. 7p holds in particular if 

dyXs = — — ; — - — dzX2 
Jc + e^m 
em 




Jc + e'^m 
eAs 



dyX2 



After some calculations we can prove that the functions 

A2 = cos a Vi{y, z, e) + sina ^(y, z,e) + a J K{x, e)dx 

A3 = — \ — i^[cosaV2{y,z,d) - sina \/i(?/, 2;, e) \K{x,e), 
Jc + "in V / 



e mx 



Jc + e^m 



are solutions of (16. 7p . where K is an arbitrary function and Vi, V2 are 
functions which satisfy the Cauchy-Riemann conditions: 

dyVi{y,z,e) =dzV2{y,z,e) 
d^Vi{y,z,e) = -dyV2{y, z, e). 



29 



Corollary 6.1 The all trajectories of the Chapliguin skate (e = 0) under 
the action of the potential field of force with components (0, mg, 0) can be 
obtained from the Descartes approach 

Proof. 

In fact, for the case when e = the classical approach for Chapliguin- 
Carathodory's sleigh gives the following equations of motion 

x = 

y = mg + sin xfx 
z = — cosxn 
sin xy — cos xz = 

Hence, by derivation we obtain 



d 



dt sinx' 



gcosx 



as a consequence. 



( x = Co, Co ^ 
gsinx 



y 

z 



' Co 
,gsmx 



+ Ci) cosx 
+ Ci) sinx 







or, 



x = 

y = {gt cosxo + Ci) cosxq 
z = {gt cos Xq + Ci) sin xq 

Clearly, the solutions of these equations coincide with the solutions of (I6.6p . (16. 7p 
with the subsidiary conditions 



are particular cases of the equations obtained from the Descartes ap- 
proach. 

Corollary 6.2 The all trajectories of Chapliguin -Carathodory's sleigh by 
inertia can be obtained from the Descartes approach 
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Proof 

Let us suppose that 



Xj = Xj{x,e), j = l,2 



then the all trajectories of the equation (16.61) can be obtained from the for- 
mula 

r. {\2{x, e) cos X — eAs sin x)dx 

y = yo + J 



Z = Zo-J 

t = to + J 



A3(x, e) 

{X2ix, y, z, e) sinx — eX^ cosx)dx 



dx 



A3(x, e) 



A.Jx, e) 



On the other hand, for the Chapliguin- Caratheodory sleigh by inertia 
from the classical approach we deduce the following equations 

JcX = €(1 

my = sinx/U 
mz = — cosxfi 
ex + sin xy — cos xz = 

Hence, after straightforward calculations we obtain the system 

in 

x = qCo cosfgex + C), = — 

Jc + me^ 

y = Co(sin(gex + C) cos x — qe cos{qex + C) sin x) 
z = Co(sin(gea: + C) sin x + qe cos(ge + C) cos x) 

m 



where q^ = , 
Jc + w^e 

which are particular case of the equations fl6.6p with 

A2 = Cq sin(gea: + C) , A3 = Coq cos{qex + C) 
Evident that in this case 



2||v| 



{Jc + me )A3(x, e ) + mX2{x, e ) = mCg 
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7 The rigid body around a fixed point in the 
Suslov and Veselov cases. 



In this section we study one classical problem of non-holonomic dynamics 
formulated by Suslov ^12j . We consider the rotational motion of a rigid body 
around a fixed point and subject to the non-holonomic constraints (a, uj) = 
where u is a. body angular velocity and a is a constant vector. Suppose 
the body rotates in an force field with potential ?7(7i, 72, 73)- Applying the 
method of Lagrange multipliers we write the equations of motion in the form 

lu = [lu X ^] + [7 X ^] + /x§ 

7 = [7 X cj] (7.1) 
(a,a;) = 

Where 

/ = diag{Ii,l2,h) 
7 = (71 = sin z sin x, 72 = sin z cos x, 73 = cos z) 

/i, I2, I3 are the inertial moment of the body. 

If we assume that the vector a = (0, 0, 1) [12], then 

IiUi = jsd-y^U - 72dy:,U 

I2LJ2 = lldy^U - 739^1 f/ 

{h - 12)^1^2 + 72C^7i U - 7i9^2 f/ + yu = 2^ 

71 = -73^^2 

72 = 73t^l 

^ 73 = ll^2 - l2^l 

The above system always has two independent first integrals 

Ki = lihuf + hujl) - [/(71, 72, 73) 
K2 = 7? + ll + li 

For the real motions K2 = 1. 

By the Jacobi's theorem about the last multiplier, if there exits a third 
independent first integral K3 which is functionally independent together with 
Ki and K2, then the Suslov problem is integrable by quadratures [12] 

To determine the integrable cases of the Suslov problem seems interesting 
the following result which we can prove after straightforward calculations. 
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Proposition 7.1 Let us suppose that the potential function U in (7.2) is 
determine as follows 



U 



1 



ihfil + hfil) - h 



where /ii, fi2(ire solutions of the partial differential equations 



73 1 



dfii dfi2 
72Tr- + 7i- 



'972 dji dj3 
then the equations (7^, fi7.3\ ) admits the first integrals 

hUi = IJ,2, h^2 = -1^1 



0, 



(7.3) 



(7.4) 



(7.5) 



The aim of this apartat is to propose the Descartes approach for heavy rigid 
body in the Suslov case. 

Let us suppose that Q = S0{3), with the Riemann metric 



G 



Is h cos z 

/s COS z (/i sin^ X + I2 cos^ x) sin^ z + L^, cos^ z (/i — I2) sin x cos x sin z 
(Ji — 12) sinxcosxsinz /i cos^ a; + /2 sin^ x 



det G = I1I2I3 sin^ z, 
In this case we have that the constraints are 

tJs = -v^ a; + cos zy = 

Hence a = (1, cosz, 0). By choosing the vector b and c as follow 

6= (0,1,0), c= (0,0,1) 

we obtain that T = 1. Consequently 

V = A2(cos zd^ - dy) - X^d^ 

The differential equations generated by v and condition (16.41) in this cases 
take the form respectively 



(7.6) 
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(a, rotv) = -vv- dzP2 — dyPs + cos zdxPs = 



(7.7) 



After the change 71 = sin 2; sin x, 72 = sin 2; cos x, 73 = cos 2; the sys- 
tem (I7.6P and condition (17.71) can be written as follow 



' 71 = 77k/^i73 



72 



hh 
h 

hh 



^273 



(7.8) 



sin 2(73 



where 



^ 972 ^71 ■ 



72Tr— + 7it/— ) — co?>xdyjji2 — suixdyjjii = (7.9) 



'^73 



973 



fi2 = — /i(cosa;A3 + sina;A2), 
Ati = hi— sin xXs + cos XA2) 

We shall study only the case when 

Hj = iJj{x,z), j = 1,2 
Hence, we obtain the equation (17. 4p . 

Corollary 7.1 The function Hi, ^2 '■ 

5^(71,72) , ^ 2 , 2 N 

/^i = g:^^ + *l(72 + 73>7l) 

5^(71,72) , ,Tr ^ 2 , 2 N 

/^2 = 7^ + *2(7l + 73.72) 

OI2 

satisfies the equation i7-4\) . 
Corollary 7.2 Let jji, H2 are such that 

55'(7i,72) 



J = 1,2 



(7.10) 
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then the potential function ( [y.g[ j and first integrals Iji7.5\ ) are respectively 

u = ^{h{^f + H^f) - h 

2I1I2 O'Ji c'72 



dS 
_dS_ 



(7.11) 



The following particular cases produces the well known integrable cases |T2] : 
The Suslov, Kharlamova-Zabelina and Kozlov subcase. 
The Suslov Subcase 
If 

S = Ci7i + C272, Cj = const, j = 1,2 

then 

1^1 = Ci, /i2 = C2 

U = const. 

which correspond to the Suslov subcase. 

The integration of the equations (17. 8p in this case produces the following 
solutions 

C2 C\ 

Ui = — , UJ2 = - — 

Cih . . , yZ/f + , /2C2COS/3 



^iici + /|C| ^ hh 

72 = , ^^^^ Sin /3 Sin ^ V^ZM + M^ + 



/1/2 



hCi cos /3 



//2(72 j2(^2 

73 = sin/?cos((W^^—- ^^t + a) 
V 

where Ci, C2, a, /3, are the arbitrary real constants. 
The Kharlamova-Zabelina Subcase 
If 



S 



/i + Ci7i + ^272)^ + 



CC2I2 



where h, Ci, C2, C are arbitrary constants, then 



71 ■ 



CCih 



Cj^/i + C2/2 



72 
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IiCi + I2C2 



h + Ci7i + C272 + 



CC2I2 



C]^/i + C2/2 



C2 



= TTTltT^ V ^ + ^^^^ + ^^^^ - c^gfc 
^ t/ = /i + Ci7i + C272 

As a consequence we deduce the Kharlamova-Zabelina subcase [6]. 
The solutions of the equation (17.21) . (17.81) are 



I 2^ 2 



C2 



h + C^l. + C2^. + ^^^^^, 



h + Ci7i + C272 



a/ /iCf + /2CI 

U = h + Ci7i + C272 

7j = aj (r - + &j(r - C4) + = 7j(^, C'l, C2, C3, C4), J = 1,2 

73 = - 7?(r, Ci, C2, C3, C4) - 7l(^, C-i, C2, C3, C4) = VA(r,Ci,C2,C3,C4) 



where 



/1/2 



-P4(''"5 Cl) C2, C3, C4) 



4 



C I1I2C1C2 



hI,C, 



hCl + -^2^*2 



2 ' 



P4 is a polynomial of four degree in r. 
Kozlov Subcase 

If we suppose that Ii = I2 and 



S = -2Carctan J + / ^(7? + 71)^(7? + li) 

{D{u)y = — 



where h and C are arbitrary real constant, hence, 

72^^ 



fJ'2 

u -- 



^2 + ^2+7i^(7r + 72') 

+ 721^(7? + 72^) 
-h + - if - ll = -/i + 73 
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which correspond to the Kozlov subcase. 
The equations fl7.8l) in this case take the form: 



X 



C cos z 
sin^ z 



sm z 



which are easy to integrate. 

The solutions of the equation of motions are: 



X = xo + C J 

y = yo-c J 



t = tn 



(1-7|)2D(1-7|) 
dl3 

(l-7imi-7|) 
dj3 



Xo + C 

yo-c 



(7.12) 



73^^73 



v/(l- 71)^4(73, /i,C) 

^73 

v/(l-7i)A(73,/^,C^) 



[ P4(73, h, C) ^ h^l - 2jI - 2h^l + 273 + /i - 



Corollary 7.3 Let fii, /i2 are the functions: 

fii = ^i(7i + 7|,7i) 
fi2 = ^2(7? + 7|,72) 

then the solutions of ^7.8^ are the following functions: 



where 



73 = Vl-7?(^) -7i(i-) 

f. dr 
t = to + j 



(r-ro), J = 1,2 



Fi(7i) = ^1(71 + 72,71)1 
^2(72) = *2(7? + 7|,72)| 
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7?+7|=l-7| 



As a particular case we obtain the Tisserand Subcase. 
Tisserand Subcase 



The interesting solution of the equation (17.41) are 



f /ii = \ /hi + ai(7| + iD + hi-fl + hill) = + 73, 7i) 

I = V/l2 + 02(7! + 7?) + b2ll + /2(72) = ^2(7? + 7I, 72) 



which produce the following potential function U : 

U = /l/ll+/2/i2 + (/l&l + /2a2)7l' + (^iai+^2&2)72 + (^l«l+^2a2)7|+A/l(7l)+^2/2(72) 

where aj, bj, hj, j = 1,2 are arbitrary real constants and fj, j = 1,2 are 
arbitrary functions. 

The case when fji'jj) = aj-jj, j = 1,2 was studied in [H], where aj, j = 
1, 2 are real constants. 

The case when = 0, j = 1, 2 is well known as Tisserands case [12]. 

After integration the equation (17. 8p in the Tisserand case we obtain the 
following solutions 



Heavy rigid body in the Veselov case 

In this example we study the problem of non-holonomic dynamics for- 
mulated by Veselov in which in certain sense is opposite to the Suslov 
problem. In this problem we consider the rotational motion of a rigid body 
around a fixed point and subject to the non-holonomic constraints 



Suppose the body rotates in an force field with potential ^7(7i, 72, 73)- 
Applying the method of Lagrange multipliers we write the equations of mo- 
tion in the form 




(7, u) = y + cos zx = 



luj = [Ilj X o;] + [7 X 
7 = [7 X ^] 



dU 
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where / is a matrix such that / = diag{Ii, I2, h)- 

The Descartes approach for this system produces the following equations: 




and 

^ _ ^ + cos z(— - — ) - 
dx dz dz dy 

where 

Pi = I3 sin^ ZX2 

P2 = {h — h + {h — h) cos^ x) cos z sin^ ZX2 + {h — h) cos x sin x sin 2; A3 
Ps = {h sin^ X + Ii cos^ x) A3 + {I2 — h) sin x cos x sin z cos 2;A2 

Finally it is interesting to observe that the construction the Descartes 
approach for the Federov case [14], i.e., 

(w,7) = a 

it is necessary in the above example make the change y ^ y + at, a = const.. 
Hence we obtain that that the equations generated by the Descartes vector 
field are 

X = X2 

y = — cos ZX2 + a 
i = A3 



and 



dp3 dpi dp2 dps 



where 



Pi = I3 sin^ ZX2 

P2 = (-^3 — h + {h — h) cos^ x) cos z sin^ ZX2 + {h — h) cos x sin x sin 2;A3+ 
a((/i sin^ X + I2 cos^ x) sin^ z + I3 cos^ z) 

P3 = {h sin^ X + Ii cos^ A3 + {I2 — h) sin x cos x sin 2; cos ZX2 

(7.14) 
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The equation 07.131) can be represented as follow 

sin^z{l3 sin^ z + cos^ z{Ii sin^ x + I2 cos^ x)dzX2{x, z) + 
{I2 sin^ X + Ii cos^ x)dxX3{x, z) + 
< cos X sin X sin 2; cos z{Ii — l2){dz\-i{x, z) — dxX2{.x, z)) + 

sin z cos -2(3(/i — I2) sin^ 2; cos^ x + 3/i — /a) sin^ z + I1 + l2)\2{x, z) 

cos X sin x(— 2 + cos^ — l2)Xz{x, z) + 2a sin z cos^ z{Ii sin^ x + /2 cos^ x) = 0. 

(7.15) 



Proposition 7.2 Lei A2, A3 are the solutions of the linear partial differential 
equations (7J^ then the functions 



Ui = 72 
UJ2 



A3 



- 7173A2 - ail 



UJ3 = sin 2;A2 — 073 



072 



(7.16) 



are the first integral of the equations of the rigid body in the Veselov-Fedorov 
case. 

In particular if 

h = I2, 

then the solutions of the above equation are 
A3 = sin zC{z) 



' sin^ z + cos^ ZI2 sin^ z\ = afi(z) + K{x) 
where C{z) and K{x) are arbitrary functions and 



n{z) = 



I2 \/ {I3 sin^ z + cos^ 2/2) (2 sin z — sin 5z + sin 3z)dz 



{—I3 cos Az + 4/3 cos 22; — 3/3 — J2 + /2 cos Az) 
Hence, from (17.161) we obtain 



^1 = 72C(2:) - 7i73- 



an{z) + K{x) 
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/3 sin^ z + cos^ ZI2 sin^ z 
aQ{z) + K(x) 
^2 = -7i<^(2;) - 7273 , , , = 072 



^3 



an{z) + K{x) 
a/ /3 sin^ z + cos^ 2;/2 



a/ J3 sin^ 2; + cos^ ZI2 sin^ 2; 



073 



(7.17) 
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Thus, we easily deduce the relation 

(/g sin^ z + cos^ zl2){uj'i + a73)2 = {aVL{z) + Kf. 

In particular if a = and K{x) = C\ = const then we obtain the well known 
first integral in the Veselov case [3] 
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